Effects of dissipation on a quantum critical point with disorder 
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We study the effects of dissipation on a disordered quantum phase transition with 0{N) order 
parameter symmetry by applying a strong-disorder renormalization group to the Landau-Ginzburg- 
Wilson field theory of the problem. We find that Ohmic dissipation results in a non-perturbative 
infinite-randomness critical point with unconventional activated dynamical scaling while superohmic 
damping leads to conventional behavior. We discuss applications to the superconductor- metal tran- 
sition in nanowires and to Hertz' theory of the itinerant antiferromagnetic transition. 
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The low-tcmperaturc properties of quantum many- 
particle systems are often sensitive to small amounts of 
impurities or defects. Close to quantum phase transitions 
(QPTs), the interplay between quantum fluctuations and 
random fluctuations due to disorder can destabilize the 
conventional critical behavior, leading to exotic phenom- 
ena such as quantum Grifflths effects [1, 2] and infinite- 
randomness critical points [3] as well as smeared phase 
transitions [4] (for a recent review see, e.g., Rcf. [5]). 

In particular, the QPTs in disordered quantum Ising 
magnets are governed by infinite-randomness critical 
points [3, 6] which display slow activated dynamical scal- 
ing. In a dissipative environment, the dynamics becomes 
even slower. In the experimentally relevant case of Ohmic 
dissipation, the tunneling of sufficiently large droplets 
(the ones normally responsible for Griffiths phenomena) 
is completely suppressed [7, 8]. As a result, the sharp 
quantum phase transition is destroyed by smearing [4] . 

In contrast, in dissipationless systems with continuous 
0{N) order parameter symmetry, disorder does not in- 
duce exotic infinite-randomness behavior in dimensions 
d > 1 [9] . This changes in the presence of Ohmic dissi- 
pation. It was recently shown that large locally ordered 
droplets are not frozen (in contrast to the Ising case. In- 
stead they display the exponentially slow dynamics asso- 
ciated with a quantum Griffiths phase [10] . This leads to 
the important question of whether the QPTs of continu- 
ous symmetry order parameters with Ohmic dissipation 
are also of infinite-randomness type. 

In this Letter, we answer this question and elucidate 
the nature of the transition by applying a strong-disorder 
renormalization group (RG) to the Landau-Ginzburg- 
Wilson (LGW) order-parameter field theory of the prob- 
lem. Our results are summarized as follows: The QPT is 
controlled by an exotic infinite-randomness fixed point in 
the universality class of the random transverse-field Ising 
model. The dynamical scaling is activated rather than 
power-law, i.e., correlation time r and correlation length 
^ are related via Inr ~ with ip the tunneling expo- 
nent. With decreasing temperature, the order param- 
eter susceptibility diverges as x ~ [ln(l/T)]2'^~'*/''^/T, 
and the specific heat vanishes as C ~ [ln(l/T)]~''/''^. 
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FIG. 1: (color online) Temperature-coupling phase diagram 
for Ohmic dissipation. IRFP denotes the infinite-randomness 
critical point. The phase boundary (solid) and the crossover 
line (dashed) between the quantum critical and quantum 
paramagnetic (QPM) regions take unusual exponential forms 
leading to a wide quantum critical region. Both phases con- 
tain Griffiths regions near the IRFP [10]. 



Here, (f> is the cluster size exponent. Close to the QPT, 
the finite-temperature phase boundary takes the unusual 
form Tc ^ exp(— const |r| "'"'') with r the dimensionless 
distance from the QPT and v the correlation length expo- 
nent. The exponents tp, 0, and i/ are universal and iden- 
tical to those of the random transverse-field Ising model. 
The resulting phase diagram is shown in Fig. 1. 

Our starting point is a quantum LGW free energy func- 
tional for an A'^-component (N > 1) order parameter ip 
in d dimensions. The clean action reads 



S 



dydx ip{x)T{x,y)ip{y) 



u 

2N 



dx^\x), (1) 



where x — (x, r) comprises imaginary time r and posi- 
tion X, j dx = / dx Jp^"^dr, and T{x,y) is the bare in- 
verse propagator (two-point vertex) whose Fourier trans- 
form reads r(q, w„) = r -I- ^^c^ + 7 \ujn^^^° with r the 
bare distance from criticality (the bare gap). ^0 is a mi- 
croscopic length scale, and ujn is a Matsubara frequency. 
The damping coefficient 7 depends on the coupling of the 
order parameter to the dissipative bath and the spectral 
density of the bath modes. We are mostly interested in 
overdamped (Ohmic) spin dynamics (zo = 2). However, 
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to demonstrate the special role of zq = 2, we also con- 
sider variable zq. Quenched disorder can be introduced 
by making the distance from criticality r a random func- 
tion of position r r + (5r(x). Analogously, disorder is 
introduced into and/or 7. 

To apply the real-space based strong-disorder RG [11, 
12], we discretize the action (1) by defining discrete co- 
ordinates Xj and rotor variables fjlr). It is important 
to note that the rotors do not describe individual micro- 
scopic degrees of freedom, they rather represent the av- 
erage order parameter in a volume AT^ large compared 
to ^0 but small compared to the correlation length ^, i.e., 
V'jW =/Aydy</'(xj +y). 

We first consider the large- limit of our action where 
all calculations can be carried out explicitly. We will later 
show that the results apply to all > 1. In the large- A^ 
limit, the discrete action reads 

-^'^'^M-'^n)J,j(t)3M , (2) 

where ji > and the nearest-neighbor interactions 
Jij > are random quantities, Eq is a microscopic en- 
ergy scale used to make the field dimensionless, and 
(t)j{ijJn) = Eq ''^ (^j(T)e"^"'^dT. The Lagrange multi- 
plier Xi enforces the large- A^ constraint ((vi'^^ (t))^) = 1 
for each order parameter component ip'^'' . The local dis- 
tance from criticality, ti ~ + , contains all single-site 
renormalizations and is thus always positive. 

The basic idea of the strong-disorder (Ma-Dasgupta- 
Hu) RG is to successively integrate out local high-energy 
degrees of freedom [3, 11, 12]. Here, the competing local 
energies are the local gaps and the interactions Jij . In 
the bare theory (2), they are independent random vari- 
ables with distributions Q(e) and P{J), respectively. 

In each RG step, we choose the largest local energy 
51 ~ max{ei, Jy}. If it is a gap, say 62, the unper- 
turbed part of the action is Sq = {T/Eo)^^^{e2 + 

72 l^^nP^'^OI^al^n)!^. The coupling of 02 to the neigh- 
boring sites, 5"! = -{T/Eo)J2j,uj„J'2]'^2{~(^n)4']i^ri), 

is treated perturbatively. Keeping only the leading 
low-frequency terms that arise in 2nd order of the cu- 
mulant expansion, we obtain new interactions S ~ 
— iT/Eo)J2uj„4>ii~'^n)Jij4'j{^n) between all sites that 
used to couple to ^2, with 

Jij = Jij H - . (3) 

At the end of the RG step, 02 is dropped from the action. 

If the largest local energy is an interaction, say J23 , we 
solve the two-site cluster 5*0 = [T / Eq)^^ X]i=2 sl^i + 



exactly while treating the interactions with all other 
sites as perturbations. The calculation is straightfor- 
ward but lengthy; details will be published elsewhere. 
For J23 ^ £2; £37 the two rotors 4>2 and 4>3 are es- 
sentially parallel; thus they can be replaced by a 
single rotor 02 with effective renormalized action 

S = (T/Eo) E^„(£2 + 72 |a;„|'/^'')|02K)p. For Ohmic 
dissipation, zq ^ 2, the renormalized gap is given by 



implying the relation 72 = 72 + 73 for the damping con- 
stants. The new rotor represents a cluster with effective 
moment (number of sites represented) 

M2 = fJ'2+ (5) 

The renormalized interactions of the new rotor with each 
of the remaining ones are given by 

J2j ~ Jij + J^j ■ (6) 

The net result of the RG step is the elimination of 
one site and the reduction of the energy scale f2 together 
with renormalizations and reconnections of the lattice. 
Since the structure of the RG recursion relations (3-6) 
is identical to that of the random transverse-field Ising 
model [3, 6, 13], we conclude that our system belongs to 
the same universality class. 

In (i = 1 , the RG step does not change the lattice topol- 
ogy. One can thus derive flow equations for the individual 
probability distributions of e and J and solve them an- 
alytically [3]. In d > 1, new couplings are generated 
in each RG step, and an analytical solution is impossi- 
ble. However, by implementing the recursion relations 
(3-6) numerically, Motrunich et al. [6] showed that there 
is a fixed point in the full joint distribution of the e and 
J that corresponds to the critical point of the system. 
In both cases, the critical point is of infinite-randomness 
type. At criticality, the distribution of the and Jij 
becomes singular and broadens without limit as fl — > 
under renormalization which also provides an a posteri- 
ori justification for using the perturbative RG recursion 
relations (3-6). One may be concerned about the initial 
stages of the RG in a weakly disordered system, where 
the strong-disorder method is not very accurate. How- 
ever, perturbative RG studies [14, 15] showed that there 
is no stable weak-disorder fixed point; instead the per- 
turbative RG shows runaway flow towards large disorder. 

We thus conclude that the infinite-randomness fixed 
point found here is universal and controls the transition 
for all nonzero disorder strength. Its critical behavior 
is characterized by three exponents "01 4>i and v. The 
tunneling exponent controls the dynamical scaling, i.e., 
the relation between length scale L and energy scale £7, 
which is of activated rather than power-law type 

ln(l/17)-L'^. (7) 
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It also controls the density uq of surviving clusters via 
nci ~ [ln(l/ri)]~'^/'^. (j) describes how the typical moment 
^ of a surviving cluster depends on Q, 

fi ^ In^ (l/n) , (8) 

while v determines how the correlation length ^ depends 
on the distance r from criticality via ^ ^ In one 

space dimension, the exponents are known exactly from 
Fisher's analytical solution [3]: 7/1 = 1/2, = (1 + ^5)72 
and 1/ = 2. In two dimensions, they were determined nu- 
merically [6, 16], yielding ip = 0.42 . . .0.6, ^ = 1.7 . . . 2.5 
and ly = 1.07 . . . 1.25. For d = 3, the scaling towards an 
infinite-randomness fixed point has been confirmed [6], 
but estimates of the exponent values are still lacking. 

We emphasize the particular role played by the Ohmic 
dissipation (zq = 2) of the magnetic modes. To this end, 
we consider how the recursion relation (4) is modified for 
zq 7^ 2. For the superohmic case, zq < 2, we find 

e~,-^=a[e^^ + er]+0{j^,^) , (9) 

where x = (2 — zq)/zq and a is a constant [17]. Thus, 
the multiplicative structure of (4) is replaced by an ad- 
ditive one. As a result, the local gaps e are only weakly 
renormalized for zq < 2. Near criticality, the distribu- 
tion of the interactions J becomes extremely singular 
while that of the gaps e remains narrow. The critical 
point is therefore not of infinite-randomness type but 
conventional with power-law scaling t ~ although 
the dynamical exponent z can become arbitrarily large 
as Zq — > 2~ . (Similar behavior was found at a percola- 
tion QPT [18].) For the subohmic case, Zq > 2, the sharp 
transition is destroyed by smearing because rare regions 
can statically order independently from each other [5, 10]. 

We now turn to the behavior of observables which is 
similar to the random transverse-field Ising model [3, 6]. 
However, there are a few differences caused by the or- 
der parameter symmetry and the damping of the modes. 
Summing over all surviving clusters using (7) and (8) 
gives unusual scaling forms for the order parameter sus- 
ceptibility X and the specific heat, 

Xir, T) = i [ln(l/r)]^*-'^/^^ e, {r'^^ ln(l/r)) , (10) 

C(r, T) = [ln(l/r)]-''/'^ Qc (r"'^ ln(l/T)) , (11) 

where and Oc are universal scaling functions. At 
criticahty, this leads to C - [ln(l/r)]~'^/'^ and x ~ 
[ln(l/T)]2'^-''/'''/r. The dynamic order parameter sus- 
ceptibility at criticality can be derived similarly. On 
the real frequency axis, '^xi^ + iO) ^ [ln(l/w)]2'^-''/'^/tj. 
This implies that low-temperature inelastic scattering ex- 
periments at the location of the order parameter Bragg 
peak should see a sharp upturn in the scattering inten- 
sity ~ [ln(l/w)]2'?^-''/'^/a; with w ^ 0. The scaling form 
(10) of the susceptibility can also be used to infer the 



shape of the phase boundary close to the QPT. The finite- 
temperature transition corresponds to a singularity in 
0^(a;) at some nonzero argument Xc- This yields the un- 
usual form Tc ~ exp(— const |'''|~'"^) shown in Fig. 1. The 
crossover line between the quantum critical and quantum 
paramagnetic regions displays analogous behavior. 

The infinite randomness at criticality leads to peculiar 
behavior of the correlation functions. The average corre- 
lation function G(x) is dominated by the rare events of 
two distant sites belonging to the same surviving cluster. 
This yields [3, 6] G (x) ~ |x|-2(<i-'^V')_ in contrast, a typ- 
ical pair of sites is not in the same cluster, and develops 
exponentially weak correlations, — InGtyp (x) ~ |x|'''. 

Our explicit calculations are for the large-iV limit of 
the 0{N) LGW theory. To discuss their relevance for 
finite N, we contrast the cases of Ising {N = 1) and con- 
tinuous (N > 1) symmetries. In the former, sufficiently 
strong Ohmic dissipation freezes the magnetic droplets 
(the localization transition in the dissipative two-state 
system [19]) leading to a destruction of the sharp transi- 
tion by smearing [4]. Recently, this was confirmed in 
a numerical strong-disorder RG [20]. In contrast, for 

> 1, isolated droplets continue to fluctuate but with a 
tunneling rate (gap) that is exponentially small in their 
size [10] because O(A^) clusters are right at the lower crit- 
ical dimension of the transition. This exponential size 
dependence of the gap requires the multiplicative struc- 
ture of the recursion (4) for the merging of two rotors. 
We conclude that this multiplicative structure is valid for 
all N > 1. Since (3) just reflects standard perturbation 
theory, all our RG recursion relations and the resulting 
infinite-randomness critical point apply to the general 
0{N) case with A^ > 1. This has been confirmed for un- 
damped dynamics {zq = 1) in the case of 0(2) symmetry 
and for a general O(A^) rotor model in Refs. [17]. 

In the remaining paragraphs, we summarize our re- 
sults, we discuss applications, and we consider open ques- 
tions. We have studied the effects of dissipation on the 
quantum phase transition in a quenched disordered sys- 
tem with 0{N) symmetric order parameter. For Ohmic 
dissipation, we have found an infinite-randomness fixed 
point while the behavior for the superohmic case (includ- 
ing undamped dynamics) is conventional. For subohmic 
dissipation, the quantum phase transition is destroyed 
by smearing. This must be contrasted with the case of 
Ising symmetry for which undamped dynamics already 
leads to an infinite randomness fixed point [3, 6] while 
Ohmic dissipation causes a smeared transition [4]. All 
of these results are in agreement with a general classi- 
fication [5] of phase transitions in the presence of weak 
disorder based on the effective dimensionality of rare re- 
gions: If finite-size regions are exactly at the lower criti- 
cal dimension, the critical point is of infinite-randomness 
type. If they are below the lower critical dimension, the 
behavior is conventional; and if they can order (freeze) 
independently, the transition is smeared. 
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Our theory has several apphcations. For instance, 
the superconductor-metal quantum phase transition ob- 
served in thin nanowires [21] was studied using a LGW 
theory analogous to (1) in one dimension with an 0(2) 
(complex) order parameter and Ohmic dissipation [22]. 
The effects of disorder on the thermodynamics of this 
problem are described by our theory. Transport prop- 
erties can also be calculated using the methods of Refs. 
[23]. Analogously, our theory should apply to arrays of 
resistively shunted Josephson junctions. 

A second potential application is the Hertz-Millis the- 
ory [24, 25] of the (incommensurate) itinerant antifer- 
romagnetic transition. In this theory, a LGW free en- 
ergy analogous to (1) is derived from a microscopic elec- 
tron Hamiltonian H = i/band + -f^int + -ffdis consisting 
of a nontrivial band structure -ffband, a Hubbard-like 
interaction Hmt and a random potential TJdis by inte- 
grating out the fermionic degrees of freedom in favor of 
the bosonic order-parameter field. The validity of the 
Hertz-Millis approach to these transitions is still a con- 
troversial question as several experiments, in particular 
in heavy fcrmion materials [26] , have shown marked dif- 
ferences from the predicted behavior. Disorder effects 
are a much-discussed possible reason for these discrep- 
ancies [26, 27]. Our theory provides explicit results on 
how dissipation and disorder can yield activated dynam- 
ics, quantum Griffiths phenomena, and non-Fermi liquid 
behavior. This should make an experimental verification 
or falsification of the disorder scenario much easier. Note 
that attention must be paid to the long-range RKKY part 
of the interaction neglected in (1). It can produce an ex- 
tra subohmic dissipation of locally ordered clusters [28] 
which leads to freezing into a "cluster glass" phase [29] at 
a low non-universal temperature Tcg determined by the 
strength of the subleading RKKY interactions. However, 
the infinite-randomness fixed point controls observables 
in the broad quantum critical region above. 
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